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We study the synchronization of a hnear array of globally coupled identical logistic 
maps. We consider a time-delayed coupling that takes into account the finite velocity 
of propagation of the interactions. We find globally synchronized states in which the 
elements of the array evolve along a periodic orbit of the uncoupled map, while 
the spatial correlation along the array is such that an individual map sees all other 
maps in his present, current, state. For values of the nonlinear parameter such that 
the uncoupled maps are chaotic, time-delayed mutual coupling suppress the chaotic 
behavior by stabilizing a periodic orbit which is unstable for the uncoupled maps. 
The stability analysis of the synchronized state allows us to calculate the range of 
the coupling strength in which global synchronization can be obtained. 
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I. INTRODUCTION 

Coupled oscillator models are widely used to model complex dynamics in non-equilibrium 
extended systems, and their synchronization has attracted a lot of attention in recent years 
Q]. In studies of coupled ensembles of nonlinear oscillators, different situations have been 
considered (identical or non-identical units, periodic or chaotic single-unit behavior, local 
or global coupling) , and a rich variety of synchronization phenomena has been found (for a 
recent review, see j2[). 

In the field of coupled map lattices, the paradigmatic model, originally introduced by 



Kaneko [3|, y| , is the ensemble of logistic maps with mean field global coupling: 

N 

x,{t + 1) = (1 - e)/[x,(t)] + ^ 5^ (1) 

i G [1, A^], f{x) = ax{l — x), and e is the coupling strength. For relatively large coupling 
global (full) synchronization occurs: the array synchronizes on the manifold Xi = ... = Xn, 
where the dynamics of an element is generated by the uncoupled map. For weaker coupling 
cluster (or partial) synchronization occurs: the array splits into K clusters of Ni,...,Nk 
elements mutually synchronized 0, Q . 

A characteristic of many biological and physical systems is time-delayed coupling in the 
interaction among many units. In the case of globally coupled units, the introduction of time 
delays makes the spatial coordinates of an element relevant in spite of the infinite range of the 
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for one-dimensional arrays 



mean-field interaction. This situation was considered in Ref. 
of coupled phase oscillators. It was shown that in the limit of short delays the ensemble 
approaches a state of frequency synchronization, and that this state might develop a spatial 
nontrivial distribution of phases. In two-dimensional arrays, distance-dependent time delays 
induce a variety of patters including traveling rolls, steady patterns, spirals, and targets jsj . 
Here we study retardation effects in a linear array of logistic maps: 

N 

N 



x,it + 1) = (1 - e)/[x.(t)] + ^Y1 /t^^-(^ - ^^^•)]' (2) 



where Tij = k\i — j\ is proportional to the distance between the ith and jth maps and k is 
the inverse of the velocity of the signal that travels through the array. In a previous work 
we considered the case in which the uncoupled maps evolve in a periodic orbit of period 2 
(when 3 < a < 1 + v^) • We found that for weak coupling the array divides into clusters, and 
the behavior of the individual elements within each cluster depend on the delay times. For 
strong enough coupling global synchronization occurs, where the dynamics of an element is 
periodic of period 2, generated by the uncoupled logistic map. The spatial correlation of 
the elements along the array is such that if k is even, at time t all elements are in the same 
state, while if k is odd, at time t neighboring elements are in different states. In both cases 
an individual map sees all other maps in his present, current, state. 

In this paper we extend the previous study and consider that the uncoupled maps can be 
either periodic or chaotic (i.e., 3 < a < 4). We find that for adequate couphng strength and 
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time delay, global synchronization occurs. In the globally synchronized state all elements 
evolve along a periodic orbit of the uncoupled logistic map. Remarkably, this orbit might 
be unstable for the uncoupled maps. In particular, when the uncoupled maps are chaotic, 
time-delayed coupling might suppress chaos, stabilizing an unstable periodic orbit. For small 
arrays we study the stability of the globally synchronized solution and calculate the minimum 
coupling strength above which the unstable orbit of the uncoupled maps becomes stable for 
the time-delayed coupled maps. The numerical simulations are in excellent agreement with 
the stability analysis. 

This paper is organized as follows. In Sec. II we analyze the existence and the stability of 
the globally synchronized state. In Sec. Ill we present results of the numerical simulations 
and the stability analysis. Finally, in Sec. IV we present a summary and the conclusions. 

II. GLOBALLY SYNCHRONIZED SOLUTIONS 

A special class of solutions of Eq. (j21) is characterized by the fact that, for all pairs i, 
j, the signal received by map i at each time corresponds to a delayed state of map j that 
coincides with the present state of map i: 

Xj{t~Tij) =Xi{t). (3) 

Thus, each element "perceives" the array as being fully synchronized, in spite of the fact 
that the simultaneous states of different elements might not coincide. In these globally 
synchronized solutions each element evolves along a limit cycle of period P of the uncoupled 
logistic map with a given phase, such that we can write 

Xi{t) = xo{t + (j)i) (4) 

with Xo(t) a particular realization of the limit cycle, used as a reference orbit. The condition 
for this solution to satisfy the evolution equation is 

(pi - (pj + rriijP = Tij = k\i - j\ (5) 

for all i and j, where rriij are arbitrary integer numbers. The symmetry of the delays, 
Tij = Tji, implies that 



0i - (pj + rriijP = (f)j - (f)i + rUjiP. 



(6) 
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Thus, the phase differences (pi — (f)j cannot be arbitrary but have to be either 0j — (f)j = UijP, 
or (pi — (j)j = P/2 + riijP, with integer number. 

We shall refer to solutions with mod (0j — (f)j,P) = V i and j as in-phase solutions, and 
solutions with mod (0j+i — (f>i,P) = -P/2 V i as anti-phase solutions. Since mod — 
(pi, P) is an integer number, the period P of the orbit for the anti-phase solution has to be 
even. The in-phase and anti-phase solutions verify Eq. (0) only for certain delay times. For 
the in-phase solution, 

mod {(pi - (p^, P) = mod {k\i -j\,P) = Q (7) 

V i and j only if /c = nP with n an integer number; for the anti-phase solution, 

mod (0,+! -(Pi,P)= mod {k, P) = P/2 (8) 

only if A; = P/2 + nP with n an integer number. 

The existence of these globally synchronized states is independent of the coupling 
strength; the only requirement is that the periodic orbit is a solution (stable or unstable) of 
the logistic map. 

To analyze the stability of the globally synchronized solutions we turn the delayed Eq. 
into a non-delayed equation by the introduction of auxiliary variables: 

Uimit) = Xi{t- m) (9) 

where 1 < i < N and < m < M with M = max(rjj). In terms of these new variables Eq. 
(2) becomes 

yi,m-i{t) si m ^ 0, 

y^m{t + l)=< ^ (10) 

[ (1 - e)/[yio(t)] + jf EjLi f[yj,k\j-i\\ si m = 0. 

Next we define the vector 

z = (z/io, 2/20, • • • , ym] yu, 2/21, • • • , ym] yiM, y2M • • • , 2/7vm) (h) 

which has 7V(M -|- 1) components. The anti-phase solutions of period 2 can be written as 

Z = (Xq,, Xb, . . . , Xf), Xa, • • • ), ~ (-^b, -^a, • • • , ^ay 3^6, • • • ) (l^) 

and the in-phase solutions of period 2 as 

Zj = [Xa, Xa, • • • ; Xh, Xh, ■ ■ ■), Zj = (Xfc, Xh, ■ ■ ■ ] X^, X^, ■ ■ ■) (1^) 
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where Xa and are the points of the period 2 orbit of the logistic map. We re-write Eq. 
(7) as 

Zi{t+l)=Fi[zi{t),...,ZMiM+l){t)\ (14) 

The in-phase and anti-phase solutions are fixed points of F^: 



F{F{Z]:\)) = F{ZfX) = Z\% 



2,1 \ ^1,2 



(15) 



To analyze the stability of these solutions we need to calculate the eigenvalues of the matrix 



OF, 



dzh 



dF, 



(16) 

This matrix has N'^{M + 1)^ elements which are either 1, 0, [1 — (A^ — l)e/N]f'{xa,b) or 
e/Nf'{xa,b)- Even for a small array size there is a large number of variables involved and 
the eigenvalues have to be calculated numerically. We observe that each matrix in the r.h.s 
of the last equation has (M + 1) x (M + 1) blcks with dimensions N x N. Denoting these 
blocks as J-'ij, with i,j = 0, ...M, it is easy to see using Eq. ^Jthat the all blocks Tij with 
i > have all entries 0, except the blocks J-'i+i^i which are the N x N identity, 1^. We also 
note that J^oo = [1 — {N — l)e/N]f'{xa,b)lN while the J-'oj are nondiagonal, for example J^oi 
is 

/ e/Nf'ixa,b) ... 
e/Nf'{xb,a) ... 








e/Nf{x,,b) 
e/Nf{xb,a) 



(17) 



III. RESULTS 



In this section we present numerical simulations and results of the stability analysis, that 
demonstrate global synchronization in the in-phase and anti-phase solutions discussed in 
the previous section. The stability analysis can only be done for small arrays and small 
delay times, since the size of the matrix A [Eq. (13)] increases as kN"^. For large arrays 
and/or large delays, we simulate the equation (0). To solve the delay equation Q we need 
to specify the evolution of Xi{t) at times 1 < t < max(rjj). We evaluated this by taking 
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a random number ranging from to 1 and by letting the array evolve initially without 
coupling. 

First, we show results for the anti-phase solution, which exists for k even. For k = 1, 
we find that for all values of a there is a value of e above which the anti-phase solution of 
period 2 is stable. Figure [T] shows the absolute value of the maximum eigenvalue, |Amaa;|, 
as a function of e, for an array of = 12 maps and three different values of the parameter 
a. For a =3.5 (dot-dashed line) the maps without coupling evolve in a limit cycle of period 
4, for a =3.83 (dashed line) the maps without coupling evolve in a limit cycle of period 3, 
and for a = 4 (solid line) the maps without coupling are chaotic. For clarity the dotted line 
indicates the stability boundary |Amax| = 1- In the three cases, for large enough coupling 
the anti-phase solution of period 2 is stable (|Amaa;| < !)• Notice that the coupling strength 
above which the solution is stable increases with a. 

We verified numerically that for larger arrays the anti-phase solution is stable. Figure 
lll^a) displays, as an example, a bifurcation diagram for = 50 (a = 4 and k = 1). The 
bifurcation diagram is done in the following way: we chose the same initial condition for all 
values of e, and we plot the 100 time-consecutive values Xi{t) (with t large enough) for a given 
element i of the array. Figure |H^b) displays the same but for a neighboring element. Above 
a certain coupling strength the array synchronizes in the period 2 orbit of the uncoupled 
map and the bifurcation diagram for the two elements coincide. The synchronization in the 
period 2 orbit is surprising since for a = 4 the period 2 orbit is unstable for the uncoupled 
maps. While the anti-phase solution is stable for e > 0.6 (see Fig. Q), Fig. |21 shows that 
the array synchronizes in this solution for a slightly larger coupling strength (e ~ 0.7). The 
critical coupling strength, ecru, above which global synchronization occurs depends slightly 
on the initial condition, and increases with a and A^. Figure 01 displays the critical value of e 
(calculated averaging over 100 different initial conditions) vs. a. ecru increases linearly in the 
parameter region where the uncoupled maps are periodic, and abruptly in the parameter 
region where the uncoupled maps are chaotic. Figure IH (solid line) shows that ecru also 
increases with the system size A^. 

Notice that below the critical coupling strength, ecru, the bifurcation diagrams shown 
in Figs. El^a) and Efb) differ. This is due to the fact that for e < ecru the array splits 
into a complex clustered structure. The clustering behavior in the simpler case when the 
uncoupled maps evolve in a period 2 orbit was studied in ^ . 
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For larger time delays and k odd, the interval of coupling strength in which the anti-phase 
solution of period 2 is stable becomes more narrow. As an example, Fig. Eldisplays |Amaa;| vs 
e for a = 3.5 and A; = 1, 3, 5, and 7. Note that in a wide range of coupling strength |Amaa;| is 
slightly larger than 1. In this parameter region, starting from random initial conditions there 
is a transient time in which the array approaches the anti-phase solution; after this transient 
the array exhibits a complex spatio-temporal behavior. The transient time increases with 
A^; as an example. Fig. El displays the mean value, < x >= J2iLi^i^ time, for four 
different system sizes. The study of this unexpected effect of the system size is the object 
of future work. 

Next, we show results for the in-phase solution, which exists for k = nP. Figure [3 
displays the bifurcation diagram for two elements of the array, and a = 3.5, k = 4, and 

= 50 (for a = 3.5 the uncoupled maps evolve in a orbit of period 4). We observe that 
above a critical coupling strength {ecru ~ 0.23) the array synchronizes in the period 4 orbit 
of the uncoupled map. As in the case of the anti-phase solution, for coupling strengths below 
ecrit the bifurcation diagrams shown in Figs. El^a) andEl^b), differ. This is due to the fact 
that the two elements belong to different clusters. The dashed line in Fig. HI shows that ecrit 
increases with the system size A^. 

For arbitrary values of k, a, and e we found a rich variety of complex spatiotemporal 
behaviors. The characterization of the different dynamic regimes is the object of future 
work. 

IV. SUMMARY AND CONCLUSIONS 

We studied the synchronization of a linear array of identical logistic maps. We consider 
time-delayed mutual coupling with delay times Tij that are proportional to the distance 
between the maps {rij = k\i — Depending on the time delays and on the coupling 
strength, different synchronization regimes might occur. If the coupling is weak the array 
usually splits into a complex clustered structure. If the coupling is large enough, global 
synchronization occurs. In the globally synchronized state each element of the array sees 
all other elements in its present state [xi{t) = Xj{t — rjj)Vz, j], and all the elements of the 
array evolve along a periodic orbit of the uncoupled maps. The spatial correlation along the 
array is either periodic or homogeneous depending on k. If k is odd the array synchronizes 
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in anti-phase, such that the state at time t of two consecutive elements is Xi{t) — Xo{t), 
Xi+i{t) — Xo{t + P/2) (where XQ{t) is a particular realization of the orbit of period P, used 
as a reference). If = nP the array synchronizes in in-phase, such that the state at time 
t is Xi{t) = Xo(t) Vz. For parameter values such that the uncoupled maps arc chaotic, 
mutual delayed coupling suppresses chaos rendering the evolution of the elements of the 
array periodic in time. Thus, an important consequence of our analysis is that delayed 
coupling might allow controlling a chaotic array by rending an unstable periodic orbit of the 
uncoupled maps, stable. 
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FIG. 1: Stability analysis of the anti-phase solution of period 2 for k = 1 and N = 12. We plot 
the largest eigenvalue of the matrix A [Eq. (13)] vs. the coupling strength for a =3.5 (dot-dashed 
line), a =3.83 (dashed line), and a = 4 (solid line). 
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FIG. 2: Bifurcation diagram obtained numerically, integrating Eq.(2) with a = 4, k = 1, and 

N = 50. We plot the values of two consecutive elements, Xi (a) and Xj+i (b). Notice that after a 
complex bifurcation scenario the two elements of the array synchronize in a period 2 orbit. 
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FIG. 3: Critical coupling strength above which synchronization occurs vs. the nonlinear parameter 
a. k = 1 and N = 100. ecrU was calculated averaging over 100 different initial conditions. 




FIG. 4: Critical coupling strength above which global synchronization occurs as a function of N 
for a = 3.5 and A; = 1 (solid line); a = 3.5 and k = 4 (dashed line). 



11 



1.4 




' 0.2 0.4 0.6 0.8 1 

e 

FIG. 5: Modulus of the largest eigenvalue of the matrix A as a function of e for the antiphase 
solution and a = 3.5, k = 1 (o), A; = 3 (x), k = b (*), and /c = 7 (+). 
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FIG. 6: Temporal evolution of the mean value < Xi{t) > for four different system sizes, N = 12 
(a), = 30 (b), TV = 50 (c) and AT = 80 (d). The parameters are a = 3.5, = 5, and e = 0.6 
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FIG. 7: Bifurcation diagram obtained numerically, integrating Eq.(2) with a = 3.5, A; = 4, and 
N = 50. We plot the values of two different elements of the array, Xj (a) and Xj (b). Notice that 
after a period-halving bifurcation the elements of the array synchronize in a period 4 orbit. 



